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Abstract. We have initiated the first CO interferometer survey of early-type spiral galaxies (SO-Sab). We observed five early- 
type spiral galaxies with HII nuclei (indicating circumnuclear starburst activities). These observations indicate gas masses for 
the central kiloparsec of ~ 1 — 5% of the dynamical masses. Such low gas mass fractions suggest that large-scale gravitational 
instability in the gas is unlikely to be the driving cause for the starburst activities. We estimated Toomre Q values and found that 
these galaxies have Q > 1 (mostly > 3) within the central kiloparsec, indicating that the gas disks are globally gravitationally 
stable. From the brightness temperatures of the CO emission we estimated the area filling factor of the gas disks within the 
central kiloparsec to be about 0.05. This small value indicates the existence of lumpy structure, i.e. molecular clouds, in the 
globally-gravitationally stable disks. The typical surface density of the molecular clouds is as high as ~ 3000 M Q pc~ 2 . In the 
light of these new observations, we reconsider the nature of the Toomre Q criterion, and conclude that the Toomre Q parameter 
from CO observations indicates neither star formation nor molecular cloud formation. This argument should be valid not only 
for the circumnuclear disks but also for any region in galactic disks. We tentatively explore an alternative model as an initiating 
mechanism of star formation. Cloud-cloud collisions might account for the active star formation. 
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1. Introduction 

Early-type spiral galaxies (SO-Sab) generally have lower 
total gas masses than late-type spiral galaxies (Sb-Scd), 
and their global star formation ra tes (SFR) are generally 
less than those of late- type spirals jYoung & Scovillel 1 1 99 it 
iRoberts & Havnesl Il994l) . Circumnuclear starbursts (HII nu- 
clei) are also less frequent in early-type spiral galaxies than 
in late-type spiral galaxies. However, when early-type spi- 
ral galaxies do host circumnuclear starbursts, their extinction- 
corrected Ha luminosities (measuring SF activity) are much 
high er on average than those o f late-type spiral galaxies 
JHo. F ilippen ko & Sargenl Il997al iKennicutlL Il998l) . Similar 
tren ds are also seen in 10//m surve ys of nearby galaxies 
JPevereuxI 1 1 98l Idiuricin et all 1 1 994 . 

High resolution CO observations are critical for under- 
standing the active nuclear SF in early-type spiral galax- 
ies. Although several high resolution surveys of molecu- 
lar gas in galactic centers have recently been com pleted 



llSakamoto et allll 999aHSofne et allEooilHelfer et allEool 
Garcia-Burill o et all l2003h . they include very few early-type 
spiral galaxies, especially with HII nuclei. This paper presents 
the first results from a CO survey of early-type spiral galaxies 
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with HII nuclei using the Nobeyama Millimeter Array (NMA), 
specifically addressing the region of the nuclear starburst activ- 
ity. 



2. Samples and CO Observations 

Our sample of early-type spiral galaxies is selected from 
the op tical-spectroscopic survey of Ho. Filipp enko & Sar gent 
( 1997a) with the criteria: (a) morphological type of SO-Sab, (b) 
hosting an HH-nucleus [Note that HII nucleus here means a cir- 
cumnuclear region with SF activity but without AGN], (c) not 
edge-on [i < 70°], (d) distance < 25 Mpc, and (e) no obvious 
interaction. From the complete sample of 1 1 galaxies satisfy- 
ing these conditions, we have observed 5 galaxies for the first 
time and found one already observed and reported (Table 
[The other five galaxies are NGC 3073, 4245, 4424, 4448, and 
4694.] 

Table [2 also includes Ha equivalent widths, and star for- 
mation rates (SFR ) within the centra l 1 kpc. The SFR is cal- 
culated as in Saka moto et alJ dl999bl) from the Ha-photom etry 
with 4"x2" aperture by Ho, Filippenko & Sargent ( 1997b). We 
made the extinction correction using the Ha/H/3 line intensity 
ratio, calculated the Ha luminosity from the central kiloparsec 
region assuming constant line intensity in the area, and used the 
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Ha-to-SFR conversion formula ( iKennicuttl Il998l) . This pro- 
cedure usually gives SFRs consisten t with those based on in- 
frared observations dSakamoto et all Tl999bh . There are, how- 
ever, some caveats for these optical SFRs. They are likely to be 
lower limits, because some Ha-emission from embedded SF 
in dust may have been missed, and because absorption in the 
background stellar continuum may lead to underestimation of 
the Ha-luminosity. The HII nuclei do not host AGNs by defi- 
nition, and therefore there is no contamination of AGNs in the 
estimation of SFR. 

Table 1 

12 CO(/ = 1 - 0)-line observations were obtained using 
NMA in November, 2002 through May, 2003. Each object was 
observed for two 10-hour tracks in the C- and D-array con- 
figurations of NMA. The complex gain was calibrated every 
20 minutes using nearby quasars. The fluxes of the quasars 
were calculated against URANUS; uncertainties in the overall 
flux scales are ~ 15%. The visibility d ata were calibra t ed us- 
ing the UVPROC-II software package ( Tsutsumi et al. I ll997l) 
and mapped with CLEAN in the NRAO/AIPS package. The 
resultant 3-D cubes have resolutions of 4 - 5" (natural weight- 
ing) and rms noise of ~ 1 7 mJy /beam in a 1 0.4 km s~ 1 channel. 
(Full channel maps and kinematic maps will be presented in our 
full survey paper). These interferometric maps include > 85% 
of the total single dish CO line fluxes (Okuda, in preparation) 
for the central 15" region in each galaxy. 

3. Molecular Gas in Starbursting Nuclear Disks 

We detect significant CO emission in the circumnuclear regions 
of all 5 early-type spiral galaxies. Figure^shows CO integrated 
intensity maps and velocity fi elds [NGC 3593, the six th galaxy 
in our sample, is reported in ISakamoto et alJ dl999al) l. In all 
6 galaxies, the CO emission is strongly peaked in the central 
region at R < 500pc. Linear scales of 1 kpc along major and 
minor axes are indicated by the crosses in Fig.n In Tabled 
we include estimates for the dynamical mass (Md yn ), gas mass 
(M gas ), and average gas surface density (£ gas ) within a radius of 
500 pc from the dynamical center of each galaxy. The gas mass 
is estimated using CO fluxes within the central 1 kpc aperture 
on the galactic disks, and a Galactic CO-to-H2 conversion fac- 
tor of X co = 3.0 x lO^cirT^Kkms-'r 1 and a factor of 1.36 
to account for He and the other elements. The dynamical mass 
is calculated using the formula Md yn = RV 2 /G. The velocity 
at R = 500 pc is estimated from the velocity difference be- 
tween R = +500 pc along the major axis in position-velocity 
diagrams. 

Figure 1 

The circumnuclear regions in these galaxies are regarded as 
starburst regions by two definitions. First, the gas consumption 
timescales, M gas /SFR ~ 2 x 10 8 yr on aver age, are as short as 
those of starburst galaxies (Kennicutt, 1998). Second, the cur- 
rent SFRs ~ 0.7M o yr _1 are larger than past-averaged SFRs, 
i. e. M star /10 10 yr ~ 0.3M G yr _1 (assuming M star ~ M dyn ). In 
addition, the Ha equivalent widths, ~ 36 A, are much larger 



than typical valu es ~ - 10 A for early-type spiral galaxies 
jKennicutdl 1 998h . indicating that their current SFRs are higher 
than most other galaxies. 

Despite the starbursts, the average gas surface densities 
at R < 500 pc are lower in the early-type spiral galaxies, 
~ 150 M pc 2 , than in late-type spir als with and wit hout HII 
nuclei (mostly > 200 M pc~ 2 , see ISakamoto et ail Il999bl) . 
Even the non-starbursting interacting galaxy NGC 5195 has a 
much higher surface density of 500 M pc~ 2 at R < 500 pc 
jKohnoet all [2002). Thus, the large-scale gas surface density 
is not the only determinant of nuclear SF. These global surface 
densities are very close to, or higher than, the average surface 
density of molecular c louds, 170 M pc 2 , in the Galactic disk 
( Sol omon et allll987l) . They, however, are not a sufficient con- 
dition for initiating SF. 

Figure |2] shows the molecular gas masses plotted against 
dynamical masses at R < 500 pc. Late-type spir al galax- 
ies (Sbc-S c d) w ith HII nuclei dSakamoto et all Gi>99b; 
ISofue et all 12003) are also plotted. The solid lines indicate 
constant gas-to-dynamical mass ratios. The early-type spiral 
galaxies are distributed between M gas /Mdy n = 0.01 and 0.1 
(mostly in the range 0.02-0.05). Evidence has accumulated for 
a lower conversio n facto r Zm i n metal rich regions such as 
galactic centers dWilsonl ll995UArimoto etall ll996). If we 
adopt an Xco that is three times as small as the assumed one, 
Mgas/Mdyn becomes three times as low. Given these low gas 
mass fractions, it is unlikely that the circumnuclear gas in these 
early-type spiral galaxies has undergone large-scale gravita- 
tion al collapse and that this initiated th e nuclear starbursts 
(see Lfog & Solomonlll984tlWadaLll992l) . Our conclusion that 
this sample of early-type spiral g alaxies has low M, as /M dvn is 
very different from the result bv ISakamoto et all dl99 9b) who 
found that the later-type spiral galaxies with HII nuclei mostly 
have Mg as /Mdyn > 0.1 (see Fig.|2j. 

In order to quantify the gravitational stability of the 
gas disks, we estimate the Toomre Q parameter averaged 
within central 1 kpc. We adopt the formula suggested in 
Sakamoto et al. ( 1999b) with a small modification. Q is de- 
scribed as 

. (M^Mdyn)'^ 

Q = a - , (1) 

M ga s 

where M lT = Rcr 2 /G and cr is the velocity dispersion. We as- 
sumed a rotation curve of V oc i.e. f3 = 1 for rigid rotation 
and /3 — for flat rotation, and then, a = J2(j3 + 1). We notice 
that a (and thus Q) is insensitive to the rotation curve, and be- 
tween 1.4 (flat rotation) and 2 (rigid rotation). The dotted lines 
in Fig. |2 show constant Q values in assumi ng cr = 10km s _1 
and a — 2 (similar to Sa kamoto et allfl999bl) . As already sug- 
gested by the small values of M gas /Mdyn, all the early-type spi- 
ral galaxies have Q > 1 (mostly Q > 3), indicating that the 
circumnuclear gas disks are globally gravitationally-stable de- 
spite their clear starburst activity. We note that cr is likely to 
be la rger than our assumed value of 10 km s _1 ( Ke nnev et all 
Il993h. If we adopt cr = 30 km s~ 1 and three times as small Xco 
dArimoto et all Il996l) . the Qs become nine times as large, re- 
sulting in Q > 27 (i.e. much more stable). 
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Figure 2 

The peak brightness temperature 7b in the channel maps is 
1 -2 K. Assuming that the CO gas is optically thick as observed 
in molecular clouds in the Galaxy dScoville & SandersH l987) 
and that the gas excitation temperature T ex is similar to the 
dust temperature of 30 K in the Galactic center dGautier et all 
the molecular gas area filling factor is / ~ 0.05 
(=7b/T ex - 7cmbX where the cosmic microwave background 
temperature is Tqmb ~ 3 K. The filling factor could be lower 
if the dust temperature is higher. Therefore, the molecular gas 
is not uniformly distributed but already confined in molecu- 
lar clouds, although the gas disk is globally gravitationally- 
stable. The typical gas surface density of the clouds becomes 
as high as ~ 3000 M Q pc~ 2 (= 1.^/f), This is much larger 
than the average value, 170 M pc~ 2 , for molecular clouds in 

the Galactic disk and close t o the value, 250 Mb pc~ 2 , for the 

J 1 1 — — I 

ones in the Galactic center (Oka et al., 2001, calculated as in 

Solomon et al. 1987). The presence of molecular clouds in a 
stable disk is understood, if the disk was gravitationally unsta- 
ble in the past and heated up after molecular cloud formation. 
Gravitational interactions among clouds may suffice to increase 
cloud-cloud velocity d ispersions and heat up the global disk 
Jjog & Ostrikeilll988h . 



4. Discussion 

Although we initially expected to detect significantly large gas 
masses so that gravitational instability within the stellar po- 
tential might account for the starbursting, our survey reveals 
that the early-type spiral galaxies with circumnuclear starbursts 
have extremely low M gas /Md yn in the central kpc. Thus large- 
scale gravitational instability is not a likely cause of the star- 
bursts. Although Toomre Q-type gravitational instability is of- 
ten used for interpreting SF in galactic disks, its relevance 
for the cii'cumnucle ar SF of early- type spiral galaxies is un- 
supported (see also, iThornlev & WilsorAll99"5l:IWong & BlitzL 
l2002l iBoissier et all l20oC These results may suggest recon- 
sidering the nature of the Toomre Q. 

If CO is detected, molecular clouds must already exist un- 
less the mass function of molecular cloud s is extreme l y dif- 

i 1 r — — 1 

ferent from th ose in the Galactic disk (Scoville et al., 1987; 

Solomon et all 11987) and in the Galactic center ( O ka et all 



2001 ); we have shown the existence of dense molecular clouds 



in our sample. Hence Q from CO observations does not indicate 
molecular cloud formation. In Q-\alue calculations, a velocity 
dispersion among molecular clouds ~ 10 km s -1 (an effective 
sound speed c s for ~ 10 4 K) is usually adopted, rather than 
the molecular gas temperature (~ 10 K, or c s ~ 0.3 km s~'), al- 
though the gas temperature should be more relevant for the col- 
lapse of protostellar cores. Thus, Q cannot indicate SF directly. 
Instead, Q obtained in earlier studies might indicate the forma- 
tion of giant molecular cloud association (GMA), i.e. a group 
of molecular clouds. Therefore, GMA formation has been im- 
plicitly assumed to be necessary and sufficient conditions for 
SF in earlier Toomre Q arguments. Considering this assump- 
tion, we are not confident that the Toomre Q-type gravitational 
instability is relevant for SF. [Note, however, that a low Q-value 



means a relatively high gas surface density. Thus, active SF is 
more likely to be found in low Q-\alue regions if it is initiated 
by some other mechanisms than the Toomre 2-type instability.] 
The above argument should be valid not only for circumnuclear 
disks but also for any region in galactic disks. Perhaps the in- 
ternal structure of molecular clouds, e.g. lumpiness and/or tur- 
bulence, is more important for star formation than large-scale 
gravitational instability. 

One would think that massive bulges in early-type spiral 
galaxies could be a key to understanding the g alaxy-type de- 
pendence of circumnuclear SF activity. In fact, de Jong ( 1996) 
found that the bulge effective radius is 100-300 pc, while the 
disk scalelength is 1-5 kpc. Hence bulge potentials could sig- 
nificantly affect gas dynamics within the central kpc, and cir- 
cumnuclear rotation curves would be flatter in early-type spiral 
galaxies than in late-type spirals. A flat rotation curve causes 
a high shear velocity due to differential rotation. For exam- 
ple, assuming an angular velocity Q = 1 kms -1 pc~', the ex- 
tension of a molecular cloud a = 20 pc and the amplitude 
of the epicyclic motion / = lOpc (« ct/k; random velocity 
cr ~ 10 km s _1 and epicyclic frequency k ~ 1 AO.), the shear ve- 
locity, v = (£2 - dV/dR){a + I), becomes as high as 30 km s in 
a flat rotation curve, but km s _I in a rigid rotation curve. Thus 
the shear velocity dominates the random velocity dispersion in 
a flat rotation curve. These facts, in addition to the existence of 
high density molecular clouds in our sample, leads us to spec- 
ulate on cloud-cloud collisions induced by galactic shear as a 
cause of starbursting in early-type spiral galaxies. Obviously, 
cloud collision will enhance lumpiness and star formation in 
molecular clouds, although the molecular clouds may not con- 
tract globally in the strong shear regions. 

In the Galaxy there is considerable evidence that the ex- 
isten ce of molecular clouds is necessary but not sufficie nt for 
SF jMoonev & Solomorl Il988l IScoville & Goodt Il989l) . and 



that cloud -cloud collisions m ight initiate active SF in many 
instances JLorerl 19761: ISc oville. Sanders & ClemensL ^1986; 



Odenwal d et all 1 1 992f lHasegaw a et all 1 1 9941) . For a thin disk 
(2-D), the timescale for cloud collision is f co n = 1 /nav, where 
n = Sg as /Mc is the typical surface number density of molecu- 
lar clouds, Mc and a are typical cloud mass and diameter, and 
v is the velocity of a test cloud relative to field clouds. We here 
adopt M c = 10 6 M Q , a = 20 pc, E gas = 150 M Q pc" 2 . These pa- 
rameters indicate that the surface density of individual molec- 
ular clouds is 3000 M pc" 2 , which is close to the observed 
value. If we assume a shear velocity of v = 30kms"', the col- 
lisional timescale becomes f co u ~ 10 7 yr. Given this timescale, 
the observed SFR of ~ 1 Moyr" 1 is achieved even if the SF 
efficiency per collision is as low as 10 %. This is consistent 
with the low SF efficiencies in molecular clouds observed in 
the Galaxy. Therefore, the cloud collision model induced by 
galactic shear can reasonably explain the observed SF prop- 
erties in our sample. Note that the assumed cloud has a tidal 
diameter of ~ 20 - 150 pc at R = 100 - 500 pc, and thus sur- 
vives against stellar gravitational potentials. This tentative idea 
will be revisited in our full survey paper. 

The results from our first CO survey of early-type spiral 
galaxies require a mechanism for triggering active SF in low 
gas density regions. The mechanism cannot be the Toomre Q- 



4 



J. Koda et al.: Starbursting Nuclear CO Disks of Early-type Spiral Galaxies 



type gravitational instability, and we have discussed a tenta- 
tive model based on cloud-cloud collisions. We are extend- 
ing the CO survey with the Nobeyama Millimeter Array and 
Nobeyama 45m telescope. We will observe all early-type spi- 
ral gala xies with central starburst in lHo. Filippenko & Sareent 
(1997a). In addition to the CO survey, higher resolution obser- 
vations of the distribution of SF regions within the central kpc 
at infrared wavelengths, and observations of physical condi- 
tions in the gas at some other molecular lines, e.g. rare CO iso- 
topes and HCN (denser gas tracer, i.e. proximate sites of SF), 
will be complementary to fully understand the circumnuclear 
SF. 
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Name Morph. 


d 

(Mpc) 


< T 
(mag) 


D 25 

(kpc) 


EW(Har) 
(A) 


SFR 500 
(M G yr-') 


(10 7 M o ) 


■ Wirt 
M, 

clyn 

(10 9 M o ) 


r500 
J gas 

(%) 


v300 

gas 

(M G pc 


~ 2 ) Reference 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 972 SAab 


21.4 


-20.17 


22.10 


34.22 


1.58 


42.6 


4.5 


9 


542 


This work 


NGC 3593 SA(s)0/a: 


5.5 


-17.25 


7.84 


13.13 


0.28 


12.9 


2.4 


5 


164 


KS99 


NGC 3729 SB(r)apec 


17 


-19.24 


13.95 


76.74 


0.99 


13.9 


2.7 


5 


177 


This work 


NGC 4064 SB(s)a:pec 


16.9 


-19.39 


21.97 


55.73 


0.43 


3.1 


1.1 


3 


39 


This work 


NGC 4274 (R)SB(r)ab 


9.7 


-19.35 


19.53 


4.65 




8.6 


5.1 


2 


109 


This work 


NGC 4369 (R)SA(rs)a 


21.6 


-19.40 


13.13 


29.76 


0.39 


10.5 


4.6 


2 


134 


This work 



Note: (1) Galaxy name; (2) Morphological type; (3) Distance based 
on Virgo infall model of iTullv & Shaval (1984); (4) Total abso- 
lute 6-band magnitude corrected for Galactic and internal red- 
dening and redshift from Ho et al. (1997b); (5) Face-on isophotal 
diameter at fi^ = 25 mag arcsec 2 in the 6-band, corrected for ex- 
tinction and inclination from RC3 (de Vaucouleurs et al. 1991); 
(6) Ho- equivalent width; (7) Star formation rate within 500 pc 
of the nucleus; (8) Gas mass within 500 pc of the nucleus, from 
CO interferometer observations assuming a conversion factor of 
X co = 3.0 x 10 20 cnr 2 (Kkms~'r'; (9) Dynamical mass within 
500 pc of the nucleus; (10) Gas-to-dynamical mass ratio within 
500 pc from the center; (11) Surface gas density averaged within 
500 pc from the center; (12) Reference for CO interferometer ob- 
servations. KS99: Sakamoto et al. (1999b). 
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Fig. 1. Upper panels: CO integrated intensity maps. Contours 
represent (0.1, 0.2, 0.3, 0.4, 0.9, 1.0) x peak inten- 
sity. The peak intensities are 58.2, 27.6, 14.9, 26.5, 8.57 
Jy beairr'km s _1 from left to right. Circles and crosses indi- 
cate the fields of view (65") and their centers in the observa- 
tions, respectively. Crosses at the lower right corners represent 
1 kpc along major and minor axes. No primary beam correc- 
tion has been applied to these maps. Lower panels: velocity 
field maps. These maps are made using > 3cr emissions in each 
10.4 km s~' channel. Contour intervals are lOkms 1 for NGC 
4369 and 20kms 1 for the others. Dark sides are receding. 
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M dyn (R<500pc) [M ] 

Fig. 2. Molecular gas mass vs. dynamical mass within 500 
pc of the nucleus. A Galactic CO-to-H2 conversion factor of 
X co = 3.0 x lO^crrT^Kkms- 1 )- 1 is adopted. Solid lines in- 
dicate constant gas-to-dynamical mass ratios. Dotted lines in- 
dicate constant Toomre g-values assuming rigid rotation curve 
and velocity dispersions of cr = 10kms~'. Q > 1 indicates 
gravitationally stable, while Q < 1 indicates gravitationally un- 
stable. The g-values become larger if we assume larger cr (i.e. 
g oc cr) and smaller Xco, which are more likely in the central 
region. 



